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Abstract The continuum one-photon annihilation at ip(2S) in e + e experiment is 
studied. Such contributions to the measured final state unr° and tt + tt~ at tp(2S) mass 
are estimated by phenomenological models. It is found that these contributions must 
be taken into account in the determination of branching ratios of ip(2S) — ► lutt° and 
ip(2S) — > 7r + 7r _ , as well as other electromagnetic decay modes. The study reaches the 
<^ • conclusion that in order for BES to obtain the correct branching ratios on these decay 

. modes, at least 10pb _1 of data below the ip(2S) peak is needed. 
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1 Introduction 



The study and understanding of the decay dynamics of the charmonium states J/if>, ip{2S) 
and possibly ■0(3770), is one of the most important topics in Quantum Chromodynamics 
(QCD). There have been lots of activities in experimental part of this study in recent 
years since BES P has collected the world largest ip(2S) sample during 1993-1994 and 
1995 running years at the ip(2S) peak. Among the studies, the reconfirmation of the 
vector-pseudoscalar (VP) decay puzzle (also called "/wr puzzle") between J ftp and ip(2S) 
decays [2j first observed by MARK II jHj and the first observation of the if>(2S) decay 
suppressions in vector-tensor (VT) modes |lj compared to naive perturbative QCD pre- 
dictions are of great interest. Since then, many theoretical efforts have been put to figure 
out the possible reasons why J/ip and ip{2S) decays to some final states, especially for 
VP and VT modes significantly deviate from "15 % rule" , while many other decay modes 
agree with the predictions [5]. 

*Supported by National Natural Science Foundation of China (19991483) and 100 Talents Program 
of CAS (U-25). 
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To extend the experimental study to other decay modes of ip(2S), to further explore 
the studied modes to higher sensitivity, and to match the huge sample BES has collected 
at the J/ip peak during 1999-2001 running years P, BES has collected 14 M ip{2S) data 
during 2001-2002 running year This world largest sample makes it possible to measure 
the ip(2S) decays to electromagnetic decay level, so that one can reach the sensitivity of 
separating the QED and QCD amplitudes in ip{2S) decays. 

As for the study of ?/>(3770), as an attempt to solve the "p7r puzzle" in ip(2S) decays, it 
has been proposed that the suppression of the ifj{2S) decays to some specific final states 
is due to the mixing between S- and D-wave charmonium states, or the ip{2S) is very 
close to DD threshold so that virtual process of ip{2S) — > DD cancels lots of if>(2S) decay 
amplitudes to those suppressed channels [Zj. Furthermore, the study of the old data may 
indicate ^(3770) has a large fraction of decays into light hadrons besides DD modes, in 
contradiction to the picture that it decays to DD dominantly |Hj. To confirm this, the 
best way is to measure the DD and total ^(3770) cross sections with high precision, so 
that the difference between the two gives the contribution of the non-DD decay rate, 
while the other possible way is to measure the exclusive non-DD decay modes. 

Now we know that J/ip and ip{2S) decay into hadrons, besides those with charmonium 
in final particles, through two interactions: the one-photon electromagnetic interaction 
and three-gluon strong interaction. The amplitudes of them, in general, may have a 
relative phase. This is also true for ^(3770) decaying into light hadrons. 

In e + e~ colliding beam experiment, J/if> and ip(2S) are produced by e + e~ annihilation, 
there is inevitable the process 

e + e~ — ► 7* — > hadrons (1) 

produced simultaneously, which is indistinguishable from the hadron events from J /if; or 
ip{2S) decays. So in e + e~ experiment, any final state may come from three processes: 
the charmonium three-gluon decays, the charmonium one-photon decays, and the virtual 
photon process without going through resonance (continuum process). The study of the 
charmonium decay dynamics is to determine the decay amplitudes of the three-gluon 
and the one-photon processes. In the case of the if>(2S), the one-photon decays has a 
comparable cross section with respected to the continuum process. Moreover, for some 
strongly suppessed modes like VP and VT, the three-gluon process may also have a 
comparable cross section. In this case, one has to consider three amplitudes and the 
relative phases between any of the two amplitudes in the analysis of the experimental 
data. Fortunately, there are some decay modes of the charmonium where the strong 
decays are forbidden, only decays through one-photon annihilation is allowed, like if>{2S) 
— > uj7i°, which violates the isospin, and if)(2S) — > 7r + 7r~, which violates the G-parity. In 
such situation, only two amplitudes and one relative phase are present, this substantially 
simplifies the study, and makes the determination of the one-photon decay amplitude 
possible, provided the amplitude of the continuum process being known. 

The electromagnetic processes, such as e + e~ — > um° and e + e~ — > 7r + 7r _ , are similar 
to the process e + e~ — > in the way that there are two Feynman diagrams : the 
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continuum one-photon diagram and the ip(2S) diagram. Taking 7r + 7r final state as an 
example, two diagrams are shown in Fig. 




The observed experimental cross sections of these exclusive channels consist of three 
parts: the ip(2S) resonance, the continuum and their interference. However, unlike the 
channel which the continuum amplitude can be calculated by QED, such terms for 
uj7i° and 7r + 7r~ are to be determined by experiment, i.e., we need to measure the form 
factors of a>7r° and 7r + 7r~. 

It is important to notice that the measured cross sections at the narrow resonances are 
sensitive to the experimental conditions, but the three parts in the cross section depend 
on different aspects of the experimental details. The measured resonance cross section 
depends on the energy spread of the e + e~ collider, but such finite energy spread does not 
affect the continuum term. On the other hand, the continuum term is sensitive to the 
invariant mass cut in event selection, while such cut hardly affects the observed resonance 
cross section under practical event selection criteria. These are to be discussed in detail 
in the following sections. 

The study of the pure electromagnetic decays of the charmonium states will shed light 
on the understanding of the charmonium decay dynamics. Some theorists, like Suzuki jH] , 
have tried to probe the interference pattern between one-photon electromagnetic process 
and three-gluon process of the charmonium decays in order to solve the "p7r puzzle". 
For such analysis, the experimental information of those decays of pure electromagnetic 
processes, are of particular importance ^U] , since they supply an estimation of the elec- 
tromagnetic part in the decay modes with strong interaction. 

In this paper, we first study the experiment dependence of the cross sections, then 
give an estimation of the cross sections of the pure electromagnetic process without going 
through charmonium resonance based on the form factors, finally the minimum required 
integrated luminosity at off-peak energy point is estimated if meaningful results are ex- 
pected for the decay modes interested. 
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2 Experimentally observed cross sections 



2.1 Resonance 

The cross section of the process e + e~ — > ip(2S) — > / (where / denotes a certain kind of 
final state) is described by the Breit-Wigner formula 

i27r • r e r f 

aBw{w) = (r-M^ + r?j^' (2) 

where is the center-of-mass energy, T e and Tf are the widths of ip{2S) decaying into 
e + e _ and /, T t and M are the total width and mass of ip(2S). Taking the initial state 
radiative correction into consideration, the cross section becomes |ITJ 

(Jr.c.{W) = J dxF(x, s ) ^ 1 _ n ^ 1 _ x ^2 (Xb w( s ( 1 ~ X ))i ( 3 ) 

where s = W 2 , x m = 1 — s'/s, y/s' is the experimentally required minimum invariant 
mass of the final state / after losing energy due to multi-photon emission; F(x, s) has 
been calculated in many references fTJ^JEI] an d II(s(l— x)) is the vacuum polarization 
factor. The radiative correction in the final states are usually not considered ^3 E3- The 
reasons are twofold. In the first place, the hadronic final system is very complicated and 
since the radiative corrections depend upon the details of how the experiment is done, it 
is difficult to give a general, model-independent prescription for them. The second reason 
is that our understanding of the hadronic problem is so crude that there is no need to 
worry about the electromagnetic corrections. In any case, if we find later on that it is 
necessary to do radiative corrections to the hadronic states for some specific problem, we 
can do the calculation then, because the initial state radiative corrections and final state 
radiative corrections can be decoupled to a large extent. 

The e + e~ colliders have finite energy spread. The energy spread function G(W, W) 
is usually a Gaussian distribution : 

1 (W-W') 2 

where A is the standard deviation of the Gaussian distribution. It varies with the beam 
energy of the collider. In case of BEPC/BES, A = 1.3 MeV at the CM. energy of the 
ip(2S) JB]- It is much wider than the ip{2S) intrinsic width of (300 ± 25) keV So the 
experimentally measured resonance cross section is the radiatively corrected Breit-Wigner 
cross section folded with the energy spread function: 

oo 

a exp (W) = J dW'a r . c .(W)G(W, W), (5) 
o 

where a r , c , is defined by Eq.(3). 
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Figure 2: Inclusive hadronic cross section of ip(2S): <Jbw for Breit-Wigner cross section; 
ov.c. the cross section with radiative correction; <j exp the measured cross section on a 
collider with A = 1.3 MeV. 



In Fig. |21 three cross sections are depicted: the Breit-Wigner cross section of Eq.(2); 
the cross section after radiative correction by Eq.(3), and the experimentally measured 
cross section by a e + e _ collider with A = 1.3 MeV. In the calculation of these cross 
sections, the following parameters of '0(25') are used [IT : M = 3.68596 GeV; T t = 300 
keV; T e = 2.19 keV *. From the three curves in Fig. El we see that the radiative correction 
reduces the height of the resonance. It also shifts the peak position to above the i[>(2S) 
nominal mass. The reduction factor p and the shift of the peak A-y/i are approximately 
expressed by [TH] 




where (3 is defined as 



and 




""Throughout this paper, we use this set of data and A = 1.3 MeV for numerical calculations. 



At the ip(2S) mass, (3 « 0.0779 and 5 « 0.06. So for ^(25), the reduction factor 
p « 0.51 and the shift of the peak A-y/s maa . ~ 9 keV. The energy spread further lowers 
down and shifts the experimentally measured ip{2S) peak. In the case of a collider with 
A = 1.3 MeV, the maximum height of the ip{2S) peak becomes 640 nb, and the position 
of the peak is shifted by 0.14MeV above the ip{2S) nominal mass. 

The channel deserves special discussion here, since it is parallel to those hadronic 

channels in ip(2S) decays which only go through electromagnetic interaction, such as um° 
and 7r + 7r~. Since this is an exclusive channel, there is interference between the continuum 
and the ip{2S) amplitudes. Such interference can be seen clearly from the scan of the 
ip(2S), see Fig. El In Fig. HI the cross sections of inclusive hadrons and fi + fi~ pairs are 
depicted for comparison. Here in the calculation of radiative correction, the upper limit 
of integration x m in Eq.(3) is taken to be (1 — Am 2 /s), which means all allowed phase 
space for multi-photon emission is integrated. The peak of curve is shifted more 

than that of the inclusive hadrons, to 0.81 MeV above the if>(2S) nominal mass. 



V~ f *Born{.S) = . (9) 



2.2 One-photon continuum 

In the total experimentally measured cross section, the one-photon continuum term has 
different features from the resonance. For the channel, this term is expressed in the 

Born order as 

4-na 2 
~3s 

After radiative correction, this term depends sensitively on the upper limit of the integra- 
tion x m in Eq.(3). At the CM. energy of the ip{2S) mass, the radiatively corrected cross 
section is 8.33 nb, when x m = 1 — 4m^/s. This value means integrating all the allowed 
phase space of multi-photon emission. If x m = 0.2, which means that the final fi + fi~ 
has the minimum invariant mass of J(l — 0.2)M^(2S), after losing energy to multi-photon 
emission, then the radiatively corrected cross section is 6.25 nb. In the actual experimen- 
tal situation, the invariant mass cut or its equivalence is usually imposed to remove the 
yU + p~ from J/ 1\) decays. For the resonance term, as long as 

*m » ^, (10) 

it is insensitive to x m . Another important feature of the continuum term is that since it 
is a smooth function of s, the finite energy spread of the e + e~ collider does not change 
the measured value. 

For the inclusive hadronic final states, the continuum term is expressed in the Born 
order as 

Ana 2 

(?rBorn(s) = ~ R(s), (11) 

3s 

where R(s) is the R- value at CM. energy y/s. Since R(s) is a slowly varying function 
of energy, so qualitatively it shares the same feature of this term for final state. 

After radiative correction, it depends sensitively on the upper limit of the integration x m 



6 




Figure 3: Contributions of three parts to the cross section: dashed 
line for QED continuum (cr c ); dotted line for resonance (<r R ); 
dash dotted line for interference^); solid line for total cross 
section(cr To '). 
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Figure 4: Cross sections in the vicinity of ip{2S) for hadron (a) 
and (b) final states. The solid line with arrow indicates the 

peak position. 
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in Eq.(3), i.e., on the experimental cuts; but the finite energy spread of the e + e collider 
hardly changes the measured value. 

3 Cross sections of two exclusive electromagnetic pro- 
cesses 

The general discussions in previous section is extended to exclusive processes. For illus- 
trative purpose, we discuss two pure electromagnetic ones, which are similar to e + e~ — > 

3.1 Resonance 

In experiments, usually the branching ratio is measured by the events number of certain 
channel (Nf) divided by the total number of the resonance events(iV^(25)): 



where 07 is the measured cross section of the specific channel and a^s) the total produc- 
tion cross section of ip(2S). With the ip(2S) parameters and e + e~ collider energy spread 
used in previous section, the total cross section of ip(2S) is 640 nb (=u^(2S))- The decay 
ip{2S) — > ujtt is reported a branching ratio of (3.6 ± 1.6 ± 0.6) x 10 -5 [Ej; the other 
decay ip(2S) — > 7r + 7r~ is reported a branching ratio of (8 ± 5) x 10 -5 by DASP but a 
later result by BES is (8.4±3.5iJ;jj) x 10" 6 jH]. If we check the original references, these 
numbers actually mean that the BES measured cross section of e + e _ — > un at i/j(2S) 
mass is 2 x 10 _2 nb while for e + e~ — > 7r + 7r~ it is 5 x 10" 3 nb. 

We notice that the observed cross section of these final states are 3 to 4 orders of 
magnitude smaller than the total inclusive hadron cross section of the continuum process 
which, according to Ref. [22], is about 15 nb. So it could be that a substantial part of the 
experimentally measured cross section comes from the continuum instead of the ip(2S) 
decay. Therefore it is essential to know the production rate of ujtt and 7r + 7r~ due to the 
continuum process in order to get the correct branching ratios of 4>(2S) decays to these 
modes. 

3.2 One photon continuum 

Here we give an estimation of the possible magnitudes of e + e~ — >• 7* — > ujtt and e + e~ — > 
7* — > 7r + 7r~ by currently available phenomenological models. These two processes are 
calculated by their form factors 




N f a f 



(12) 



4ira 2 \p\ 3 




(13) 



3 ' s 3 / 2 
8vra 2 \p\ 3 



-^7T7T (^) I 1 



(14) 



3 ' s 5 /2 
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where p is the momentum of the final particles, F UJW (s) and F nn (s) are the form factors 
for un and 7r + 7r~ respectively. Phenomenological model [TU] gives 

m 2 M 2 

F^{s) = ^(°) (m 2_ s )(M2 - s y ( 15 ) 

with 

F W7r (0) = 2.3 GeV- 1 . (16) 

At E cm = 3.686 GeV, it gives 3.6 x 10~ 3 nb at Born order. It is almost unchanged by 
radiative correction if x m = 0.2 in Eq.(3)t. 

For 7r + 7r~, the vector dominance model (VDM) [2B] gives : 



Fl DM (S) = (17) 
nip — s 



At if}(2S) mass, it gives 3.2 x 10 3 nb at Born order, which is almost unchanged after 
radiative correction if x m = 0.2*. 

From above estimations, we can see that possibly some of the observed un and a 
large fraction of the observed tc + tt~ events at ip{2S) peak may come from the continuum 
process instead of the i/j(2S) decay. It should be noted that different models for the form 
factors give very different results at the ip{2S) mass region j2H|, here we only want to get 
a feeling about the magnitude of the continuum contribution. 

For -0(3770), the resonance cross section is only around 8 nb; whereas the hadron cross 
section from the continuum is about 13 nb. To determine its branching ratio of exclusive 
light hadrons, we must know the continuum cross section of these final states. 



3.3 Energy dependence of measurement 

From the foregoing discussion, we know that the behavior of two terms of the cross section 
is rather distinctive with the energy variation: the resonance changes manifestly but the 
continuum part almost imperceptibly. The relative proportion of two terms of the cross 
section, which plays an important role in determination of branching ratio, could be rather 
different at different energy. 

In actual experiment, the beam energy could drift away from the position of the 
maximum inclusive hadron cross section during a considerable long time running, which 
in turn leads to considerable variation of resonance cross section. Taking / u + fi~ final state 
as an example, if the energy drifts upwards or downwards from the ip(2S) peak position 
of the inclusive hadron by 0.5 MeV, the fi + fi~ resonance cross section changes from 4.47 
to 4.20 nb, or equivalently the variation is about 6%. 

^The radiatively corrected cross section is 3.9 x 10~ 3 nb if x m = 0.3 and 3.3 x 10~ 3 if x m = 0.1. 
*The radiatively corrected cross section is 3.4 x 10 _3 nb if x m = 0.3 and 3.0 x 10 _3 nb if x m = 0.1. 
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Furthermore, for the process whose continuum cross section is comparable with that 
of resonance, its interference cross section accounts for a considerable part in the total 
cross section, and varies violently in the vicinity of the peak position. Also for final 
state, its interference cross section changes from —0.14 nb to —0.96 or 0.66 nb respectively 
for the energy drifted by 0.5 MeV lower or higher than the hadron peak position. 

Combining the variation of the resonance and the interference, the maximum change 
of the resonance cross section can reach 1.1 nb, or 25% for channel. In other words, 

if we use the result from fi + fi~ final state as an approximation, 25% uncertainty should 
be taken into account in the branching ratio estimation for the 0.5 MeV energy drift*. 

4 Other ip(2S) decay modes 

In previous sections, we have discussed the pure electromagnetic processes in ip{2S) de- 
cays. Most other ip{2S) decays which strong interaction plays the leading role also have 
continuum amplitude contribution if the ip{2S) is produced in e + e~ collision. Although 
the amplitude maybe smaller, but for those suppressed channels, like VP and VT modes, 
the electromagnetic decays and the continuum contribution could play an important role. 
For example, now the upper limit of ip{2S) — > pn branching ratio has been pushed down 
to 2.9 x 10 -5 |2]. Under BES condition, this means that the upper limit of the cross 
section is 2 x 10~ 2 nb. The coupling of this channel to virtual photon is one third of ujtt . 
Using the estimate of previous sections, its continuum cross section could be 1.2 x 10~ 3 nb 
at '0(25'). If we push down further the upper limit of this decay, we need to take into 
account the contribution from the continuum. This is also true for the upper limit of 
ip(2S) -> K + K*(892)- + c.c. Another VP decay mode, i/j(2S) -> K*°K° + c.c, has been 
measured to have a branching ratio of (0.81 ± 0.24 ± 0.16) x 10~ 4 [21]. This means that 
BES measured cross section of this channel is 5 x 10 _2 nb at ip(2S). Theoretically, the 
coupling of this channel to virtual photon is two thirds of um°, so the estimated continuum 
cross section by the form factors of Eq. (j!5)l is 3.2 x 10~ 3 nb. This has to be considered 
in high precision experiments of the coming generation accelerator and spectrometer, like 
CLEO-c and BES-III. 

In order to know whether the observed suppression of VP and VT modes in ip(2S) de- 
cays are due to the absence of strong interaction amplitude, or the destructive interference 
between the electromagnetic and the strong amplitudes, or just an incidental destructive 
interference between these two and the continuum process in our particular experiment, 
we need to know their coupling to virtual photon. 

5 Estimation of needed data taking at BES 

As discussed in previous sections, in order to extract the branching ratios of ip{2S) to wk° 
and 7r + 7r~ from experimental data, we need the coupling of these channels to a virtual 

*If the energy drift is about 0.2 MeV, the uncertainty is only 9% correspondingly. 
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photon. This can only be done by the direct measurement of their production cross 
sections nearby but off 0(25) resonance. The same is true for -0(3770) decaying into light 
hadrons. To satisfy the need of both -0(25) and -0(3770) physics, the best energy to take 
the data should be below -0(25) peak, for example, at E cm = 3.67 GeV. Here two factors 
should be considered: one is the data-taking point should not be too faraway from the 
peak position, so that the cross section of the continuum process is almost equal to that 
at the peak position; another factor is the data-taking point could not be too close to the 
peak position, so that the resonance cross section is small enough to be neglected. Taking 
the inclusive hadron final state as an example, at E cm = 3.67 GeV, the cross section 
variation of the continuum process 

x °h(E cm = 3.67 GeV) - a h (E cm = 3.686 GeV) 

^ = a h (E cm = 3.67 GeV) * L05% - 

At the same time, the ratio of the resonance cross section to that of the continuum process 

R / resonance \ _ ( Oe+e-^{2S)^hadron \ ~ 2 18% 

V continuum /£ cm = 3.67 GeV \ <?e+e-^r ^hadron )e = 3 67 GeV ' °" 

(19) 

That is to say, at E cm = 3.67 GeV, the variation of the continuum process cross section 
and the proportion of the resonance cross section are both fairly small. 

Another reason to take the data below -0(25) peak is that -0(25) and -0(3770) are too 
close to find a suitable data-taking point between them. On one hand, -0(3770) is a wide 
resonance, the resonance effect can extend far beyond the peak position. On the other 
hand, the effect from the radiative tail of -0(25) can reach the vicinity of the -0(3770) 
peak. For example, at one r(-0(377O)) below -0(3770) peak, the resonance cross section is 
0.8 nb, which is mostly DD; the radiative tail of -0(2,5) is more than 3 nb; while the cross 
section for continuum process is about 13 nb. That is to say, the radiative tail of -0(25) 
is too large to be neglected. 

Finally, we discuss quantitatively how much data are needed. 

• The typical electromagnetic process, like -0(25) — > ujtt°, is measured to have a 
branching ratio of 3 x 10~ 5 , since the cross section by BEPC/BES at ■0(25') peak 
is 640 nb, that means the ujtt cross section is 2 x 10~ 2 nb. So for every pb -1 , we 
obtain 20 produced events. 

If we use the form factor in Eq.(13) to estimate the cross section, at 3.67GeV it is 
3.8 x 10~ 3 nb with x m = 0.2 in radiative correction. If so, for every pb -1 , we obtain 
3.8 produced events. 

• Another pure electromagnetic process ■0(25') — > 7r + 7r~ is measured to have a branch- 
ing ratio of 8 x 10~ 6 , corresponding cross section is about 5 x 10~ 3 nb. For every 
pb -1 , we obtain 5 produced events. 

If we use the form factor in Eq.(14) to estimate the cross section, at 3.67GeV it is 
3.3 x 10~ 3 nb with x m = 0.2 in radiative correction. If so, for every pb -1 , we obtain 
3.3 produced events. 
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• e + e~ — > pir is one of the most interesting processes. Each of the three final states, 
p°7r°, p + 7r~, p~7T + gets 1/9 contribution as um° from pure one-photon process. So 
for every pb -1 , we expect 6.7 produced events (which is one of the three final states). 
For continuum process, for every pb -1 , we obtain 1.3 produced events. 

From the above estimates, the integrated luminosity of the off resonance data sample 
need to be the same order of magnitude as on ip(2S) resonance. 

According to the experience of BES data taking, the typical integrated luminosity 
at ip{2S) mass is about 8 nb -1 per RUN; average 25 RUNs can be taken every day, or 
equivalently 0.2 pb -1 per day. As an estimation, using these experiment data and form 
factors, the time needed for events numbers taken at the continuum of different processes 
can be worked out and are given in Table H 



Table 1: Events number and time needed 



Events number 


1 


4 


10 


25 


50 


data taking 




1.3 


5.3 


13 


33 


66 


duration 


7T + 7T~ 


1.5 


6.1 


15 


38 


76 


(day) 


pix 


3.8 


15.4 


38 


96 


192 



In order to obtain reliable signal, and taking both statistical error and data collection 
time into consideration, we recommend that at least 10 pb -1 data at off ip(2S) peak to 
be taken. From the above estimates, it gets the statistic to match the new sample about 
20 pb -1 of ip{2S) taken by BES-II. It will serve the purpose to take into account the 
contribution from one-photon continuum in data analysis. 



6 Summary 

In this paper, we discussed extensively the properties of the observed cross section at res- 
onance in e + e~ experiment. We studied the resonance and the electromagnetic processes. 
We also pointed out the possible uncertainty due to actual measurement circumstance. 
The electromagnetic final states of um° and 7r + 7r~ are from two processes: one is ip{2S) 
decay and the other the continuum contribution. Using form factors, we found the contin- 
uum contribution is equivalent to the previously reported cross sections of the resonance 
decay. This indicates that the actual branching ratios of these electromagnetic decays 
might be different from the present reported values. 

In order to obtain the correct results of ip(2S) and ^(3770) decay, we suggest that at 
least 10 pb -1 data be taken at the continuum region [E cm = 3.67 GeV). 
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